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ABSTRACT 

A pseudo-spin model is proposed, as a means to describe some transport prop- 
erties (resistivity and Hall mobility) in Bi2Sr2{CazPri-z)Cu20s+y. Our model 
is based in a double-well potential where tunneling in a given site and interaction 
between different lattice sites are allowed only through the excited states. Doping 
of the pure system by the addition of Pr increases the ratio between the activation 
energy and the tunneling constant. The model Hamiltonian displays some features 
which are present in the hydrogen-bonded ferroelectrics. Its dynamics is treated 
in the random phase approximation and the characteristic frequency (time) is used 
in a Drude formula in order to obtain some transport properties of the system, 
namely the electric resistivity and the Hall mobility. The quantities calculated in 
this work are compared with the experimental data of B. Beschoten, S. Sadewasser, 
G. Giintherodt and C. Quitmann [Phys. Rev. Lett. 77, 1837(1996)]. 
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1. Introduction 

The controversy behind the mechanism which governs high temperature superconductivity 
(HTS) in cuprates is a subject of current interest, as can be verified in the recent literature [1,2,3]. 

Anderson [2] attributes the novel phenomenology present in cuprates materials to a second kind 
of metallic state, namely, the Luttinger liquid. Zhang [3] has proposed a SO (5) theory of cuprate 
superconductivity which considers that the phenomenology of these materials might be fundamen- 
tally due to a conflict between different kinds of order. Chianchi, Moretti and Piazza [4] have 
presented a phcnomcnological model for the formation of Cooper pairs in cuprate superconduc- 
tors: in their description, coupling is mediated by large anharmonic oscillations along c direction 
of the apical ions. On the other hand, Beschoten et al [5] have investigated the charge transport 
in insulating and superconducting samples of Bi2Sr2{CazPri^z)Cu20s+y, in terms of the Hall 
mobility. The main motivation behind Beschoten et al work [5] was that the interplay of localiza- 
tion and superconductivity is of fundamental interest in the physics of electrons (holes) in strongly 
correlated and disordered systems such as the high-temperature superconductors (HTSC). 
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The aim of this paper is to propose a phenomenological model to explain the resistiv- 
ity and the Hall mobility behaviors inferred from Bcschoten et al data [5], performed in 
Bi2Sr2{CazPri-z)Cu20s+y The strategy of our work is first to propose a model Hamiltonian 
to describe the Bi2Sr2{CazPr\-z)Cu20s,+y and then to study its dynamics in the random phase 
approximation (RPA). After doing that, we use the obtained characteristic frequency (time) in the 
Drude formula, in order to obtain the resistivity of the system. The Hall mobility will also be 
studied by using the above results. 

Gaona and Silva [6], inspired in the Gorter-Casimir (GC) two-fluid model of superconductivity 
[7], as well in the formalism used to treat the phase transitions in hydrogen-bonded ferroelectrics 
[8], have proposed an effective transverse Ising model (TIM) to describe both metallic and HTS 
superconductivity. We think that an extension of this previous model could be useful to describe 
some basic features of Beschoten et al [5] experiment. 

Let us first make some considerations about two well-established phenomenological theories 
of the superconductivity. One of the main characteristics of the superconducting state is the 
macroscopic phase coherence of the wave function. In 1950, Ginzburg and Landau [9] introduced 
their phenomenological theory of superconductivity, where the order parameter must be identified 
with the macroscopic wave function ^ . They started from the idea that 5* represents some effective 
wave function of the superconducting electrons. \E' can be normalized in such a way that is 
equal to the concentration of superconducting electrons. It seems that the pioneering theory 
to treat the thermodynamic behavior of the superconducting state was proposed by Gorter and 
Casimir (GC) [7] . The two-fluid model of Gorter and Casimir was proposed to provide a basis for 
understanding the thermodynamic of the superconducting state. In this phenomenological theory 
it is assumed that the electrons in a superconductor are divided into two interpenetrating gases: 
the superconducting electrons, which have zero entropy, and the normal electrons, which have the 
usual properties of electrons in a normal metal. The normal electrons give up the condensation 
energy (3 when they become superconducting. The authors proposed that below Tc a fraction y, 
where y depends on the temperaure T of the conduction electrons occupy a set of lower " condensed" 
energy states which are associated with the superconducting properties, while the fraction {1 — y) 
remains uncondensed. The free energy function for the whole electronic assembly is made up of 
two terms: one, JjVj is similar to the electronic free energy of a normal metal and the other, Fs, 
corresponds to the condensation energy associate with the condensed state , 



where 7 is the normal electronic specific heat coefficient, and /3 is a constant representing the 
condensation energy. Gorter and Casimir assumed that the phases were not independent, and 
proposed the following free energy function: 



Close to Tc we may express (2) in the standard Ginzburg-Landau form by expanding F (T, y)as 
a power series in y, assuming y = j^fj , and retaining the first nonvanishing terms. The model has 
acquired support by its success to interpret all properties of superconductors. In particular, the 
temperature dependence of the penetration depth is correctly predicted [10,11]. 

An inspection of the GC free energy relation (2) reveals that it contains two contributions: the 
first one is related to the property of free electrons in metals to transpose the Fermi barrier, being 
responsible for the transport properties and consistent with the Pauli exclusion principle. This 
term also contains a multiplicative factor which takes into account the fact that, below Tc, the 
number of normal electrons in the superconducting sample is a temperature dependent quantity. 
The second term of the GC free energy is responsible for the pairing of wave functions of the 



Fn = -^lT\ Fs = -P 



(1) 



F{T,y) = {l-y)^/^FN+yFs . 



(2) 
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electron condensate. This term is also temperature dependent (indeed, it is proportional to the 
number of the electrons in the condensate) and above Tc it docs not contribute to the GC free 
energy. So, the GC free energy clearly displays the competition between these two terms. 

On the other hand, a very successftil model to treat various cooperative phenomena is the Ising 
model. However, the Ising model suffers from a deficiency, namely, the impossibility to exhibit a 
proper dynamics. Perhaps the simplest extension of the Ising model with a proper dynamics is the 
transverse Ising model (TIM). The TIM has been used to describe phase transition in ferroelectrics, 
ferromagncts and cooperative Jahn Teller systems [12]. This model Hamiltonian was first proposed 
by De Gennes [13] based on arguments developed by Blinc [14] to represent the basic features of 
hydrogen-bonded ferroelectrics of the KH2PO4 family. In these systems the Ising term corresponds 
to the interaction between the protons at different lattice sites and the transverse field accoimts 
for the possibility of protons occupying one of the two minima of a double potential well in a given 
site. 

2. The Model Hamiltonian 

Taking into account the preceding considerations, we propose the following effective Hamiltonian 
as a means to describe some transport properties of the Bi2Sr2{CazPri-z)Cu20s+y We write: 

H = -A\Y.a-nY: '-=^X, - \ Y: L,^-^Z^-flZ, . (3) 

i i ij 

The main ideas that we have in mind are: 

i- ) The third term of (3) is an Ising-like term, where the operator Z is related to the wave 

function of the condensate electrons (holes), and the coupling Lij favors the pairing of electrons 
(holes) at two different lattice sites, contributing to the coherence of their wave functions below 
the critical temperature Tc- 

ii- ) In the second term of (3) the operator X is related to the wave function of the normal electrons 
(holes), where fl represents the transverse field or tunneling constant. This term, which gives 
dynamics to the model, could be thought as the kinetic contribution to the model Hamiltonian. It 
could also represent a mimic for the motion of "free" electrons (holes) through the Fermi barrier, 
in close analogy with the tunneling of protons through the double-well potential barrier in the 
KH2P04-like (hydrogen-bonded ferroelectrics) case (see [8]). However, while in KH2PO4 the proton 
tunneling occurs between the right and left sides of the potential well, in the present model we 
consider the possibility of the electron (hole) to tunnel between the two discrete values of the 
phase ip of its wave function, namely between = and ip = n. According to this point of view, 
the spontaneous symmetry breaking will occurs when the system makes the choice between one of 
these two particular phases of the wave function. 

iii- ) In the first term of (3), A represents the metal-insulator activation energy, and the operator 
C will assume the eigenvalues -1-1, when the electron (hole) occupies the ground state (insulator), 
and -1, when it occupies the excited state (metal or superconductor). We observe that the electron 
(hole) could change the phase of its wave function only in the -1 cigcnstatc of C. The interaction 
of electrons (holes) between different lattice sites, which leads to the cooperative effect of the 
superconductivity, occurs only in this last case. 

With respect to the operators which appear in the Hamiltonian (3), we would like to make the 
following remarks: X and Z will be treated as linear independent operators in the space of spin- 
1/2, described by the Pauli matrices. The operator C, which has eigenvalues -|-1 and -1, could be 
treated as a " z component " of a spin-1/2 operator in another space. In this way, the Hamiltonian 
(3) is defined in the space product of the operators C and ( X,Y,Z ). 
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Before finishing this section, it is worth to mention that a Hamiltonian of the kind described 

by (3) was used by Ohtomi and Nakano [15] to explain the role of the pressure in determining the 
order of the phase transition in hydrogen-bonded ferroelectrics. They however treated (3) only in 
the static case, namely in the mean field approximation (MFA). 



3. The Dynamics of the System 



The Dynamics of the system will be treated by applying the RPA to the Heisemberg equations 
of motion. We have (with U= 1) 



dt 



(4) 



where the operators Pi = Xi, Yi, Zi must satisfy the Pauli commutation relations, 

[Xi,Yi]=iSijZi, ... . (5) 

To solve (4) in the RPA, we write 



(6) 



where in (6) we make the replacement of {Pi)^ by a constant part (Pi), which is just the mean field 
approximation (MFA) expected value, plus a small time-dependent deviation 6 (Pi) e**^*. 

Using (5) and (6), keeping only terms which are linear in the deviation S {Pi) and by putting 
{¥) = 0, we have the following set of linearized equations: 

iivS {X) - Lo (^—tP-) (Z) 6{Y)=0 



i(jS {¥) - 



Q-Lo{^-^]{X) 



6 (Z) + Lo 



1-(C) 



{Z)S{X) = 



(7) 



The solutions of (7) are: 



and 



^2,3 — ^ 



iuj6 {Z) + nS {¥) = 



wi = 



n-Lo{^-^]{x) 



1-{C) 



which are constrained by the zero-order solution 



n-Lo{^-^]{x) 



{Z)=0. 



(8) 



(9) 



(10) 



In the following we are interested only in the characteristic frequency of the electrons (hole) in 
the non-superconducting phase, i.e., in the phase given by (Z) = . 

Putting {Z) = in 3 given by (9) and dropping the subscripts for sake of simplicity we have: 
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n-Lo{^-^]{x) 



(11) 



which leads to: 



n-Lo{ 1 (x> 



1/2 



In (12) the values of (C) and {X), in the MFA, are given by 



and 



(C) = tgh 



A-n{x) 

2k0T 



(12) 



(13) 



(14) 



4. Transport Properties 

1. Evaluation of the Resistivity 

To evaluate the resistivity we will work with the Drude formula, 

m 



(15) 



where m is the effective mass of the carriers, p their concentration, e the electronic charge and r 
the scattering time. 

Let us now make the link between the dynamics of our model and the Drude formula. Since the 
carriers concentration nnist bo an increasing function of the Ca concentration, / (z), and must be 
a maximum when all the electrons (holes) occupies the excited state ((C) = —1) and a minimum 
for the ground state ((C) = +1), we propose: 

P = nf{z)(^-^^\ (16) 
where n is the density of carriers which are able to form cooper pairs. We also assume [16] that 

where ^ab and are respectively the in-plane and the perpendicular coherence lengths. The 
number 2 corresponds to the assumption that in the superconducting phase we will have only a 
cooper pair by "unit cell" of volume ^a(,^c- 
By putting 

(18) 



2|a;±| 
in (15) we get 
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p = 



2m 



|w±l 



(19) 



and putting (12) and (16) into (19) we get: 



/8m\ { 



1/2 



nf{z){l-{c)y 



(20) 



One immediate consequence of relation (20) is that the onset of the superconductivity is reached 
when the term in the square-root goes to zero. This corresponds to the softening of the mode 
which frequency is given by equation (12). 



2. Hall Mobility 

To make a connection between our model and the Hall mobility pi, we will consider the one-band 
model [17] where the Hall effect is inverse proportional to the carrier concentration p: 



Rh = — . 

pe 

Using this assumption we will have for the mobility 

Rh (T, z) t 



11 = 



p{T,z) 2m\uj±\ 



(21) 



(22) 



where we have used (19) to obtain the right side of (22). 
Putting (12) into (22) we finally obtain 



\e\h 
2m 



1-{C) 



{X) 



-1/2 



(23) 



3. A Possible Mode of Vibration 
We also can evaluate the resistivity with an equivalent way of writing the Drude formula 

TUVf /„.s 

P=^- 24 
pe''l 

where vp is the Fermi velocity and I the mean free path of the electrons (holes). 

As a means to determine the electron (hole) mean free path we start with the classical formula 



1=^. (25) 



Wo observe that in (25) we have considered the number of scatter centers by unit of volume to 
be equal to the density of Cooper pairs given by (17). We can justify this by the argument that 
the electric conduction always occurs in a regime of charge neutrality. 

The electron (hole) mean free path, I, remains to be linked to the model's dynamics. Alterna- 
tively, we can work with the collision cross-section 7rr^ taking into account equation (25). To do 
this let us write the "resonance" condition: 
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^Kr'^ = 2n\u)±\ . (26) 

In (26), we have considered the equahty in energy between a classical harmonic oscillator of 
elastic constant K and the energy separation of a two-level quantum system, which frequencies 
ui± are given by (12). 

Initially we are going to determine the elastic constant K through the following arguments. It 
seems natural to consider that the Debye frequency lvd and the mass M of the oxygen ion must 
play an important role in the scattering of the charge carriers. So, let us write the relation: 

K = aMtol, . (27) 

We can suppose that in (27) K is an effective elastic constant, where a , an adjustable perhaps 
small pure number, could be interpreted as a mimic for a weak coupling regime of the BCS [18] 
model. 

Inserting the informations contained in the equations (25) and (26) into (24), we obtain 

/Ann^yu^^ (28) 



and putting (12) and (16) into (28) we get: 



1/2 

(29) 



f{z){l-{C)f 

Since (19) and (28) are two equivalent ways of writing the resistivity, we obtain 

K = 2TTnnvF . (30) 

On the other hand an estimate of the Fermi velocity vp was given in reference [19], in terms of 
the electron (hole) mass and the in-plane coherence ^ab, namely, 

VF = l-^. (31) 
2mU 

Putting (17), (27) and (31) into (30), we get 

^ = -i^-n) = ^. (32) 

The fact that the right side of (32) depends only on the in-plane and perpendicular coherence 
lengths and on the effective mass m of the carriers leads to the conjecture that it is possible to 
find a mechanism accounting for the elastic constant K. In the following we are going to look for 
this possibility. 

In a previous work by one of the present authors [19], an effective potential was proposed to 
describe the interaction between electron pairs (holes) relating the size of this super carrier to the 
energy of the bound state. By considering the particular situation that the particle occupies the 
minimum of this effective potential and that the motion of the pair is confined in a plane [20], 
we found that this problem could be mapped into that of the quantum mechanical behavior of 
a particle described by polar coordinates in a plane. The solution of this quantum mechanical 
problem leads to an eingevalue spectrum which is doubly degenerate, except for the ground state. 
In the first excited state we have two eingevalues of the momentum q, namely: 
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g± = ±7^ . (33) 

t,ab 

Wc may consider that a pair of electrons (holes) circling in opposite directions establishes a ring 
of charge and that an oxygen ion could experiment the electric field of these pair of carriers. This 
leads to an oscillatory motion of the oxygen ion perpendicular to the plane, and in the case of 
small oscillations, this motion could be described by a harmonic oscillator which frequency coo is 
given by: 

Mu^l = • (34) 

Now, if we make the identification of the left side of (34) with the elastic constant K given by 
(32), we obtain: 

«.= (i22)|5!. ,35) 

So, according to our model the perpendicular coherence length is inversely proportional to 
the effective mass m of the charge carriers. 



5. Comparison with the Experiments and Numerical Estimates 

Let us now compare the experimental results of Beschoten et al [5] with the present model. We 
may consider that the role played by the addition of Pr in the Bi2Sr2{CazPr\-z)Cu20%+y is to 
promote certain kind of dilution which is responsible by the decreasing of the critical temperature 
of the material. We assume that the dilution affects both the values of the " tunneling constant" 
f2, as well as the activation energy A. In order to understand this, let us look at figure 1. In this 
figure we have a double- well potential where the minimum corresponds to one of the two values of 
the discrete phases of the local wave function, = or = tt. The parameter A corresponds to 
the separation in energy between the ground and the excited states. In figure 2 we show that the 
effect of the tunneling is to make the splitting of the excited state in two sub-levels. The system 
could become superconductor as far as ri/2 > A. Therefore the critical concentration of Ca ions 
{zc), i.e., the concentration below which we have no more the superconducting state, will occur at 

n = 2A. 

Looking at the resisivity formula given by (20), we can observe that its temperature dependence 
is controlled by the parameters of the model. So for Q. > 2A, the onset of the superconductivity is 
reached when: 



2n 



(^)t=t. = . ^ , (36) 



where {X) and (C) are given by (13) and (14), respectively. 

In the numerical estimates of the resistivity (equation (20)), we have three adjustables parame- 
ters: the tunneling constant Cl, the activation energy A and the ratio between the coupling and the 
tunneling constants Lo/fl. Table 1 shows the values of the parameters adjusted to reproduce the 
experimental data findings of Beschoten et al [5] . In figure 3 we plot in a semi-logarithmic scale 
the resistivity (equation (20)) of the Bi2Sr2{CazPr\-z)Cu-20%+y as a function of the temperature 
for various values of the dilution z. We used, in evaluate n given by (17), ^ab = ISA and = 3A 
[21], and m is the electron rest mass. 
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z 



A 



0.3 
0.4 
0.5 
0.6 
0.7 
0.9 



12.62 
9.68 
7.61 
5.35 
1.31 
0.9 



24.0 
19.0 
15.0 
11.0 
2.8 
2.1 



4.2 
4.0 
3.8 
3.2 
4.1 
5.0 



Table 1: Adjusted parameters values for various Ca concentrations {£) in 
Bi2Sr2{CazPri-z)Cu20s+y (A and Cl are in units of 10~^^ J). 



A linear relationship between z and A/J7 can be infered from the data of table 1. Prom this re- 
lation we can estimate a value for the critical concentration, Zc, approximately 0.51, corresponding 
to A/n = 1/2. 

In the numerical estimates of the resistivity for the Bi2Sr2{CazPri^z)Cti20%+y, in the version 
of equation (29), we have used a as an adjustable parameter. M is the mass of the oxygen ion, 
is the Debye frequency given by kpOo I'h, where we used for the Debye temperature = SOO-ff 
[22]. We have adjusted a = 0.005, consistent with the value of Mwg obtained in (34). 

We also obtain for using equation (35), the value of approximately 2. 5 A. 

Beschoten et al [5] have investigated the inverse of the Hall mobility, in 
Bi2Sr2(CazPr\-z)Cu20%j^y cupratc. They foimd that the mobility of the conduction holes evolves 
from the insulator [z < Zc = 0.52) to the high temperature superconductor [z > Zc) by using a 
one-band model with a single scattering time in order to fit the experimental data. Their main 
result is that, even for high Ca content, z > Zc , where the system becomes a superconductor, they 
still clearly observe a minimum in the inverse mobility, although its resistivity is purely metallic 
(this indicates that, even in these samples, the superconducting wave function is formed by states 
which are spatially localized). 

However in our results we can not observe those minima in the inverse mobility (equation (23)) 
for z > Zc, although they are present for z < Zc- 

In figure 4 wc use the zeros of the inverse Hall mobility and the minimum of it ( for z < Zc) as 
a mean to construct a diagram separating the various regimes: insulating, superconducting and 
metallic. This diagram agree with the experiments (figure 4 of reference [5]) unless for the fact 
that it does not reproduce the coexistence of superconducting and localization (insulating). 



It seems that the model proposed in this paper can capture the main features which are present 
in the transport properties of the Bi2Sr2{CazPri-z)Cu20s+y In particular, we have evaluated 
the resistivity which agree with the experimental data of Beschoten et al [5] . 

From the evaluated inverse Hall mobility we extracted a diagram T (temperature) versus z 
(concentration of Ca ions) showing the various regimes of the system, namely metallic, localized 
(insulating) and superconducting regions. 



6. Concluding remarks 
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By comparing two alternative ways of writing the Drude formula of resistivity, and by finding 

a possible mode of vibration (equation (34)) very relevant to the problem, we can led to a rela- 
tion which links the effective mass of the carrier (hole) to the perpendicular correlation length 
(equation (35)). 

The present research has been supported partly by CNPq-Brazil. 
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FIG. 1. Figure 1: The double-well potential, showing the ground and first excited states. At the 
minimal, the local phase of the wave function is = or <p = tt. 
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FIG. 2. Figure 2: Splitting of excited state of the potential well, due to the tunneling term Q. When 
n/2 > A the system is in the superconductor regime (for z > Zc). 
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FIG. 3. Figure 3: Logarithm of the resistivity as a function of the temperature, for various values of 
of Ca concentration (z). Theoretical results of this work (solid lines) are compared with the experimental 
ones (circles) of Beschoten et al [5] . 
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FIG. 4. Figure 4: Inverse mobility minima corresponding temperature Tmin (open squares) and super- 
conducting transition temperature Tc (filled squares) for Bi2Sr2{CazPri-z)Cu20s+y- 
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